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1. Introduction 



The Baxter-Bazhanov-Stroganov (BBS) model [5, 6] (also called the r'-^-'-model, see e.g. 
[7, 8]) is associated to the cyclic L-operators [9, 6, 10] which act in a two-dimensional 
auxiliary space 

L.(X)= ( 1 + ^ ^ ^7^ (^i - ^i) 

^ V ~ ^ + ^^Z^^' 

where A is the spectral parameter, and at each site j = 1, . . . ,n we have five parameters 
aj, bj, Cj, dj and Xj. At each site there is also an ultra- local Weyl algebra with elements 
Uj and Vj, obeying 

UjUk = UkUj, VjVfc = VfcVj, UjVfc = LJ^^'>'VkUj, id = e^''*/^, = vf = 1, N e Z>2 . 
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Since c<j is a root of unity, the Weyl operators can be represented naturally by matrices 
acting in the tensor product (C^)*^" [1, 2]. The monodromy matrix of the model is 
defined by 

and the transfer matrix is its trace in the auxiliary space t„(A) = tr T„(A) and gives 
rise to a set of commuting non-local and non-hermitian Hamiltonians of the model: 

t„(A) = A„(A) + L'„(A) = Ho + HiA + • • • + H„_iA"-^ + H„A^ 

Commutativity follows from the intertwining of the Lk{X) by the asymmetric 6- vertex 
model i?-matrix at root of unity. 

In the present paper our focus will be on the case N = 2. As it has been shown in 
[11], in this case the BBS-model is related to a generalized Ising model with plaquette 
Boltzmann weights 

iy((7i,(72,(73,(74) = ttg (l + Z)l<i<j<4 + ^4 C^i (72(73(74^ , (3) 

subject to the free-fermion condition 04 = 012^34 — ^13024 + 014023- 

Generalizing Sklyanin's method of Separation of Variables (SoV) [12, 13, 14], in 
[1, 2] we have worked out a method to find common eigenvectors of the Hamiltonians 
H^. This proceeds in two stepsij 

• Since [B{\), B(iJ,)] — 0, the off-diagonal element B(X) of the monodromy matrix 
(2) gives rise to an auxiliary set of commuting operators h^: B{X) = hiA -|- 
h2A^ + ■ ■ • -|- h„A". We construct their common right eigenvectors I^E'p) (the left 
eigenvectors (\E'p| are obtained analogously) by an inductive procedure over the 
chain size n, starting from the one site model. The eigenvalues of -B(A) form a 
polynomial in A of degree n, and from the intertwining relations we can show that 

n-l 

B(X)\^p) = Xrou;~''°l[{X + rku;-'"^)\^p), p = (p^, pi, . . . , (4) 

k=l 

where the amplitudes (ro, ri, . . . , r„_i) can be expressed in terms of the parameters 
of the model aj, . . . , >Cj, and we can use the phases p, & I^n, oi the zeros 

= -TkUj''''', k = 1, . . . ,n - 1, (5) 

of the eigenvalue polynomial for labelling the eigenvectors. 

• Having solved the auxiliary problem, after a Fourier transformation of po to p 
{p labels the Zjy-charge sectors), the eigenvalue problem of t(A) = A{X) + D{X) 
is reduced to the solution of Baxter equations. Using that A(z/fe) and D{vk) are 
raising and lowering operators on the auxiliary states |^'p), we find that the periodic 
eigenstates |$p,e) in 

t(A)|$p,E) = t^^\X\E)\^p^^), E = {Ei,...,E„_i} 

X We consider a fixed chain length n, mostly suppressing the subscript n, which had been written 
explicitly in our previous papers, so = rn,k and pn = pn,k of [1, 2], etc. 
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with the eigenvalue polynomial t(^^(A|E) = Eq + EiX^ h-E^A" (where and 

En are directly known, see (20)), are obtained via the kernels in 

|^p,e) = 5]a;-''-^°QVlP,E) p' = (pi, . . . , p„_i). (6) 

The crucial fact is now (SoV) that after splitting off a known function f{p'), the 
n — 1- variable function Q^{p'\p, E) factorizes into single variable functions Q^{pk) 
(the pk are the components of p', we often skip the charge index p of the etc.): 

q^(pV,e) = /(pO mz}Q^(p.), 

and the Q^{pk) are determined by the Baxter equations {k — 1, . . . ,n — 1): 
t^f\v^ |E) g^(pfe) = A+(i.fe) g^(pfe + 1) + ^-^{uv^) Q^{p^ - 1). (7) 
The corresponding Baxter equations for the left periodic eigenvector read 
t^^\vkm Qlip,) = u;'^-'A-{iy,)Ql{p, + l) +uj'-''A+{ujiy,)Ql{p,-l). (8) 
The functions are defined by 

n-l 

A+(A) = {oo^/xk) (A/a;)^-'^ J] ^M/^) , ^kW = Xk (A/a;)'^"^ F„(A/c.) , (9) 

m=l 

-^m(A) = ( &m + ujam X^A) ( A + d^/ ) . (10) 

We shall not need the expression for Xk: see (43) of [2], since this will cancel in our 
final formulas. The existence of a non-trivial solution to (7), (8) is provided by a set of 
functional equations, which determines the still unknown values E. 

In [2] we calculated the action of u„, the Weyl operator at site n , on an eigenvector 
l^p) of -B(A) to have the form 

u„|*p) =^|*p) + ES^fe|*P+*) with p+'== (po,...,Pfe + l,...,Pn-i), (11) 

and g and gk are certain functions depending on the parameters aj,bj,Cj,dj, >Cj and 
the components of p. Since in [1, 2] we also found a factorized expression for the 
norm (^p|^p/), we have the framework for calculating normalized matrix elements 
(^p|u„|^p/) / (^o| ^o), where = (0, 0, . . . , 0). For calculating matrix elements between 
periodic states ($p|u„|$p/), in addition we also need the solutions Q^^\pk) of the Baxter 
equations. These are available for N = 2, and our main goal is to obtain such periodic 
matrix elements in a factorized form. We achieve this by explicitly performing the sums 
over the intermediate Z2 variables. 

This paper is organized as follows: In the following Section 2 we recall the N — 2 
spin matrix element calculated in [2] and transform it into a much more compact form by 
performing the summation over an intermediate variable, still keeping the model general 
and inhomogencous. In order to proceed beyond this result, in Section 3 we specialize the 
parameters of the BBS-model such that we get the homogenous Ising model. We discuss 
the structure of the eigenvalues and the solutions to the Baxter equations. The vanishing 
of some transfer matrix eigenvalues at the zeros of -B(A) requires us to distinguish four 
cases when solving the Baxter equations. Then in Section 4 we continue the evaluation of 
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the matrix elements of the spin operator, until we are finally able to perform the multiple 
sum over intermediate spins. The derivation of the basic formula for the multiple spin 
summation is delegated to the Appendix. Section 5 deals with the calculation of squares 
of the matrix elements. In Section 6 we include a special case excluded in the earlier 
Sections and give the final formulas for the matrix element of spin operators in terms 
of the zeros of the transfer matrix and excited quasi-momenta. Then we are ready to 
compare our result in Section 7 to a conjectured formula of A. Bugrij and O. Lisovyy 
[3, 4]. In Section 8 we apply the formulas of Section 6 to obtain the matrix element of 
(7^ for the finite quantum Ising chain in a transverse magnetic field. In Section 9 we 
give our conclusions. 



2. Spin operator matrix element for the N — 2 inhomogenous BBS model 

In [2] we derived a formula for the normalized matrix element of the spin operator 
between arbitrary states of the periodic inhomogenous N — 2 BBS-model. For N — 2 
there are two Z2-charge sectors p = 0, 1 in (6). Since u„ is anticommuting with the 
charge operator V = V1V2 ■ ■ ■ v„, its only non- vanishing matrix elements are between 
periodic states with different charge p. Let Q^^''^ and Q^^^^ be solutions to the 
Baxter equations (7), (8) and the zeros of the operator polynomial B{X), see (4). 
Then our result in (58) of [2] can be written: 



j:m{p') U(p') (^(-ir'- "'""-^;^""-'''- )+i:w)) . (12) 





u 


■n\ 


^1) 


( ^0,0 1 


^0,0 ) 



p' \ ^ " ^ k=l 

n-1 



wp') = (-!)"'■ n ,n„^'l?:,y„, ■ MP') = n '°'(«)Qr<"(«), m 

l<m ^ ' V / ;_1 



, n—1 

Ojqn, On Cn 



Rkip') = - ^ Ql^'\p, + 1) Qf'\pk)llQi''\Pi) Q^iPi) X 

X (l-^)^^p^, (14) 

with Vk — —i^k{—^y'', f — roTi- ■ - Tn-i and p' = ^^Z^ Pk- The different terms in (12) 
have the following origin: H{p') is a normalization factor due to the convenient (since 
it avoids further factors) choice of normalizing by (^0,0 1 ^0,0) of the auxiliary system. 
Here \^ p,p'J = \'^o,p'J + (-l)''|^i,p;), where \'^po,p'J is an eigenvector defined in (4), 
see also (35) in [2]. For the terms in (12) recall (11). The terms at Ro{p') correspond 
to g and go in (11). The sum over k and the expression for Rk{p') arise from the shift 
in the index pk and the coefficients g^. 

In the remaining part of this Section we now show that the sum over k can be 
performed, leading to the much simpler expression (21), (22). 
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We start rewriting the factors of the Rk{p'): 

anKCn A \ _ . (0)" 



ro x„ z/fc (z/fc + Cn) 1 2 V / 2 



2 ro x„ i/fe (i., + Cn) \^ g^(°) (p, + 1) Qf'^ (p,) 

where we define Cfe = ^'fe/(ofeXfe) and use Ajf^^(A) from (9), pointing out the exphcit 
dependence on p. For obtaining the first two hnes of (15) we used 

and A['^-'^(— z/fc) = — A^^ '*'''(— z/,fc). To get the third line of (15) we used the Baxter 
equations (7), (8), where for iV = 2 we have uo = —1 and + 1 = — 1 mod Z2: 

gf)(p.) Ar(-.,) + Af-(.,) Qf)(p, + 1) Ar(-.,) + Ar(.,) 
Ql^\pk + l) i-^r-'m^k) ' Qf)(p,) tW(-^.) 

Then, using (15), (14) becomes: 

Rkip') = ^ , ^,,lr , . ^O(P') + 

2ro X„ i/fe (i/fe + Cn) \\s+k\^h -^s) \ 

n-1 \ 

+ {-ir-H('\-u,) Qi^'\pk + 1) Ql^'\pk + 1) llQi^'^ipOQf^'^ip^n- (16) 

Ij^k / 

Now the sum over k in (12) can be performed using an identity (c.f. the Appendix of 
[2]), vahd for any polynomial f{x) of degree less than n + 1 and for any n + 1 non- 
coincident points Xk- Consider a polynomial f{x) = fn^"' + ■ ■ ■ + fo, its interpolation 
through the points xi, . . . , Xn+i, and focussing attention on the coefficient of x": 

n+l n+1 n+1 ^/ n 

fe=l s^fe k=l i-L^ky^k-Xs) 

For calculating the sum over k of the first term in the parentheses of (16), in (17) we 
take f{x) = t(°)(x) and {xi, . . ., x^-i, Xn+i) = {1^1, ■ ■ ., i^n-i, 0, -(«)■ Thus we get 

^ ^k i^k + Cn) n.^fe(^fc - ^s) " Cn n::i' (-^.) -Cn n::i'(-Cn - ^.) ' ^ ^ 

where En^ is the leading coefficient of (A) . 

For the second term of the last line of (16) we will not perform the summation over 
k directly. Instead, for each p' for which we make the summation of the first term, we 
take for the summation over k the second term of (16) corresponding to p'"^*^ (which 
entails i/^ — > —Vk) in (12). Collecting all such terms and taking into account the changes 





u 


■n 




( *o,o 1 


*o,o) 
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which come hova.H{p'^^) /H{p') we perform the summation over k. Together with (18), 
we get 

_bn_ ^ I^^(O) _ t^'m , ti'\-Cn) 

Cn Us=l i-'^s) Cn Us=l (Cn " J^s) J 

The leading and the constant coefficients of t^'^\X) can be read off directly from (1),(2): 

n n n 

i=i 1=1 1=1 

Inserting these results into (12) we find that the E^^ and E^f^ terms in (19) just 
cancel the terms of the bracket at -Ro(p') in (12), and we get simply 

^ J2 N{p')Ro{p')R{p') (21) 

^ith ^ ^ tW(-Cn) t^^HCn) .22) 

nr="/(-Cn + (-i)^'n) nr="/(Cn+(-i)^'n)- ^ ^ 

We shall write R{p') = R^^\p') + R''^\p') when we have to refer to the 
separate terms on the right hand side of (22). Despite the simple appearance, for 
the general inhomogenous N = 2 BBS-model, performing the multiple sums over the 
Z2 variables seems to be a presently impossible task. However, restricting ourselves to 
the homogenous model with the parameters satisfying 

Qj = Cj = a, bj = dj = b, Xj = 1 for j =l,...,n — l. (23) 
we are able to evaluate (21) with (13), (22) completely, as will be shown in Section 4. 

3. Homogeneous Ising model 

In all following Sections we consider only the N — 2 case of the model defined by (1) 
with (23). For a fixed chain length n, we are left with only the two parameters a, b and 
the spectral parameter A. For N = 2 we have u = —1 and we represent the Wcyl op- 
erators Ujfc, Vfc by Pauli matrices acting at the k-th site. So, now our model is defined by 

/ 1+Xaf \al{a-bal)\ 
Lk{X) = . (24) 

Fixing the spectral parameter to the value X = b/a, the L-operator (24) degenerates 
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and the transfer matrix can be put into the standard Ising form 

n n 

t{h/a) = tr Li(6/a) 1.2(6/0) • --L^ib/a) = JJ(1 + ■ h/a) ■ JJ(1 + a^.^a^ ■ ab) 

k=l k=l 

- exp (ELi K at) exp (ELi a^) , (25) 

if we use periodic boundary conditions 0"^+^ = cx^ and identify 

e-^^y =tanhK* = b/a, tanhK^ = ab . (26) 

So at A = b/a we call the model (24) the Ising model. If we don't fix the spectral 
parameter to this special value, we shall talk of the "generalized Ising model" . However, 
transfer matrix eigenstates are independent of the choice of A. 

3.1. Structure of the eigenvalues 

In [1] the eigenvalues of the transfer matrix t(A) = tr Li{\) ■■■ Ln{\) with Lk{X) 
given by (1) for N = 2 and homogeneous parameters, have been calculated from the 
functional relations. From Z2-invariancc t(A) commutes with the Z2-chargc operator 
V = erf erf ... (J^ . Since = 1, the space of eigenstates of t(A) decomposes into two 
sectors according to the eigenvalues (—1)'^' (where p = 0, 1) of V. The sector p — 
is called the NS-sector, p — 1 the R-sector. The 2" eigenvalues can be written (we 
specialize assuming (23)): 



where the quasi-momentum q in each sector takes n values q = ^ m with m integer 
(half-integer) for the R (NS)-sectors. If ciq = the quasi-momentum q is called 
unexcited, for erq = 1 it is called excited. In the NS (R) sector, the eigenstates of 
t(A) have an even (odd) number of excitations: Yiq = ("!)''• 

For q = (this occurs only for the R-sector) and for q = tt we define 

- 1 + 1 

So = ^ r , = ^— - . 28) 

The quasi-momentum q = tt is in the R sector for n even, it is in the NS sector for n 
odd. The different presence of factors (A ± Sq) and (A ± s^) in (27) for n even or odd 
often makes it necessary to consider the cases of even n and odd n separately. 
Sometimes we shall use the notation Aq :— (— l)'^''Sq. 

3.2. State vectors from Baxter equations 

In order to obtain the eigenvectors of the transfer matrix t(A), we have to solve Baxter's 
equations (7) and (8). As input we use the corresponding eigenvalues t^(A) which are 
specified by the values ciq for all q in the sector p, see (27). Solving Baxter's equations, 
we should use the values t''(±rji.) of these polynomials at the values ±rfc of the roots of 
the eigenvalue polynomials of the operator Bn{X) (which is the off-diagonal element of 
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the monodromy matrix) given by the formula (A7) of [1]. For our special parameters 
(23) and N — 2 the are simply related with the s^: 

n = Sqfc , Qk = irk/n, A; = 1, . . . , n - 1. (29) 

This means that for our special choice of parameters (23), the zeros of t'^{X) may coincide 
with the Tfe, giving rise to the vanishing of the left hand sides of (7) and (8). At the 
parameters (23) all are equal: F^(A) = F{X) and from (10) we obtain 

F(A) = 6^ - A^ xl r?^""'^ = (-l)"+*^+^F"-2(r,) . (30) 

Let us compare two sets: the set {qk} which parameterizes the roots Vk, and the set 
of all possible quasi-momenta {q}. The latter set is divided into two sub-sets: the NS 
and the R sectors. The NS-sector contains pairs of quasi-momenta {qk, —qk} for odd k 
and the R-sector includes the pairs {qk, —qk} for even k. The quasi-momentum q = 
always belongs to the R-sector, and q — n belongs to the R-sector for even n and to the 
NS-sector for odd n. 

The solutions of Baxter's equations for the case of the Ising model were found in 
[2]. Here we recall the final result. For a fixed sector p and the eigenvalue polynomial 
^(A) we have to solve n — 1 systems of Baxter's equations (7) (or (8)) numerated by the 
integers k = 1, . . . ,n — 1. With respect to these data we have to distinguish four cases: 

{-ir = {-ir: 

(i) tP{rk) ^ and tf{-rk) ^ 0: 

Qk Qk (1)- 2xkrr'Firk)- 

The other three cases occur for 

(ii) ^{rk) 7^ 0, ^{—rk) = 0: ^''(A) contains a factor (A -\- ruY (both q = ±qk not 
excited), we may normalize 

g^'» = i, Ql^{i) = 0. 

(iii) tP{rk) = 0, tP{—rk) 7^ 0: ^^{X) contains a factor (A — r^)^ (both q = ±qk are 
excited), we cannot choose Qf,' (0) = 1, but we may normalize 

Ql'» = 0, g^^(i) = 1. 

(iv) t^(rfc) = tP{—rk) = 0: t^{X) contains (A^ — r^) (either q = +qk or q = —qk 
is excited): A L'Hopital procedure as described in [2] is required (in order to obtain 
eigenvectors of the translation operator), leading to 



where 

t'W - 4W (A + (-l)'^^^s,J(A-(-l)'^^^._,J, ^(q) = a^-2a2cosq + l. (31) 
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In the next Sections we shall restrict ourselves to calculate only transitions between 
eigenvectors allowing the normalization Q^'^(O) = 1, postponing to Section 6 the 
consideration of eigenvectors for which i''(A) contains factors (A — Vkf', i.e. the 
eigenvectors involving case (iii) above. 

As already observed at the beginning of Section 2, the non-vanishing spin matrix 
elements have left and right eigenstates from different sectors. Let and t^^^ be the 
corresponding eigenvalue-polynomials. With respect to these two polynomials we define 
k e P^''^ if t'' has a factor (A + TkY, i.e. we have case (ii), 
k e V'^P^ if has a factor (A — r^Y, case (iii), and 
k E D^''^ if f has a factor (A^ ~ "^D^ i-^- have case (iv). 
By D —\T)\we denote the number of elements inD — T)^^^ U D^^\ similarly for etc. 

4. Calculation of the matrix element of in the homogeneous Ising model 

We now start to evaluate (21) with (13) and (22) in our simplified model where 

C = h/a, rl = {a^ -h''){a^'' -l)/{a^ -I). (32) 

We have to observe that in the derivation of (12) given in [2], generic BBS-parameters 
leading to tP{rj?) ^ were assumed, and the solutions to the Baxter equation were 
normalized to Qfc'^(O) = 1- As we have seen in Section 3.2, in the case (iii) this 
normalization is not possible for the special parameters (23). In order not to complicate 
the derivation, in the following part of this Section we shall simply exclude state vectors 
containing k & V, adding the changes necessary for A; e "D in Section 6. Also in this 
section we shall omit the superscripts L and R in the notations of Q^'^^''\pk) supposing 
that the left/right eigenvectors are from NS/R-sectors as they appear in (21). 

Consider Ro{p'). Always one of the factors in Qf\pi)Qf'\pi) is from case (i) above, 
and excluding / G T)^ the other is from either (ii) or (iv). So always Qf\id)Qf\id) = 1. 
For / G P we have Qf^\l)Q[^\l) = since cither Q\^\i) = or (5p^(l) = depending 
on the parity of /. So, in (21) the summation reduces to the summation over pi for 
I e V, with pi = fixed for I e T). 

4.I. Calculation of R{p') 

Let us show that a common factor can be extracted from the two terms of (22). We 
first consider the case of odd n where q = appears in the R-sector and q = tt in the 
NS-sector. 

We start with the first term i?(°)(p') in (22). Now from (27) the NS eigenvalue 
polynomial t^^\X) for odd n is 

NS, nodd: ^('^^(A) = (a^- + l)(A + (-l)--s,) [] {X + r^f [] (-^^ " ' (33) 

kef>(o) levio) 

(for even n omit the bracket with s^^) since in the NS-sector only odd k appear, and 
these fall into one of the classes (ii) and (iv), class (iii) being momentarily excluded. We 
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insert C) from (33) and decompose the denominator product over / in its even-Z 

and odd-/ parts. We write the odd part as / G P'-^' U P*-^^ since for p = in cases (ii) 
and (iv) / must be odd (recall, we still exclude case (iii)): 

^ ^aeven(-C + (-l)«n) n..^(0)(-C + r.) a,^(0)(-C+(-l)«n) 

n.„cvcn(C+(-l)^-r„) 

In the last line we put a factor unity, written as quotient of upstairs a product over 
m e T>^^^ U P^-^^ and downstairs over meven. In the T) terms we omitted the factor 
(—1)'''= since from (ii) this contributes only if — 0. Now several cancellations take 
place, resulting in 

Observe that now the p'-dependence appears only in the last product over V. This 
happens because all Z-odd terms cancel and because allows only pi = 0. In the 
denominator we use ( = b/a, (27), (29) and Yli even ^ili) = ('^^" ~ l)/!"^^ ~ 1) ^o obtain 

n,eve„(-C^ + rf) = ae - l){a'b' - l))^"-'^^' (a' - l)/(a2" - 1) . (36) 

The second term in (22) can be evaluated analogously. We insert t^^\C) from 

R, nodd: iW(A) = (a^- - 1)(A + (-l)'^°so) J] + H " (^7) 

and use -^{Qi) — {(^'^^ + + 1) to get finally for n odd: 

i?('^°^^)(p')= ((-ir(a'+i)(-c Yiievii-^rri +0 

-(-l)-(a^-l)(C +(-l)-«o) Yiievii-^Tri - 0) R (38) 
with (using also (32)) 

R = ro^«-^(«/3)-('^-^)/2a"-in.ei^ ((-1)' C + r.) , (39) 
where cx = a' (3 = l-a%\ and (-l)l»""l = (-1)<^-, (-l)l^'''l = -(-1)'^° in the 
case of odd n. 

The case of n even is less symmetric between and R^^^ since now both q = 
and q = TT are in the R sector, none of them in NS. So the term containing s^r appears 
in tW(C) instead of in t(°)(-C). Also, (-l)l^'°^l = 1, (-l)l^^'^l = -(-i)<^o+a. ^he case 
of even n. In the following products, for I odd there are n/2 values and for I even we 
have n/2 — 1 values: 

i?(neven)(^,) = (n^e^((-l)''' + C) - {-ir^^^Ulevii-^r ri - 

X (a^ - 1) (C + i-irso) (C + (-l)'^^sO aV(a/3)) R , (40) 
R = rl a-\af3Y-^/'a--' J] ^ + " (^1) 
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4.2. Calculation of U{p) ■ Ro{p') 

In this subsection we shall show that the product ■ -Ro(p') can be put into the very 
simple form (49). Let us start evaluating Ro{p'). 

At the beginning of this Section we already discussed that, if we exclude case (iii), 
I V, then g[°^(0)g[^^(0) = 1, and from (ii) if / e P we have Qf^\l)Q\^\l) = 0. 
So we have to consider only / G V for which we get from (i) and (iv) (momentarily we 
suppose that / is odd, but the result (45) is the same for even I): 



Qf{i)Q\'\i)= -(-1)'^-+' 



,(- l)"2ising,4°)(-n) (-1)"-! t(i)(-r,) 



nyl(%)F-Hn) ^ ^ ^ ' 

where in the last step we used (30). The polynomial t^^^{X) is t^^\X) given by (33) 
with the factor — rf omitted, sec (31). In the first line the factor —(^—lyii+i takes 
care whether qi or —qi is excited, the minus sign comes because Qf^^ is a left eigenvector 
component. Now since we should not use / G ^ (if present, such a term leads to a 
vanishing contribution in (21)), we have 

An 1 

-^^.n/?-3n(-^). («) 

with Si := {a'^-l){-ri + {-iy^s^){-ri + {-iy°so). Inserting (43) into (42) and using, 
recall (31), 

4n _ 1 2 

rf -rl = 2 (cosg, - cosg,) F{ri)/A{q,) , H^Zj A{q,) = = . (44) 

and the trigonometric identity 

2"-i siv?qi Yik^iicosqi - cosg^) = n(-l)'+i 

we obtain 

«!»>(!) = ^--f^ n (^) . («) 

valid both for qi in R and for qi in NS. 

Let us rewrite the ratio 5*/ /(2i sing/ F{ri)) in a convenient way. By straightforward 
use of the definitions (28), (30) and (29) we find 

^2ising,F(n) + (-l)-<'a±,, ' ' - ^ ^ ^ ' 

^2ising,F(n) n + (-l)-/3±,, ' - e'^ ^ ^ ' ^^^^ 

leading to (written such that it is valid for both p/ = and pi = 1): 
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Qf\p,) Qhp.) = (-i)<"-"" tii^ . n ti)^ , (47) 



where 

^1 = 1"^- ^ for (-1)-" = ±(-1)-- ; qi = . (48) 

I A = (-1)"° Pgt K ) , m K ) K ) 

Multiplying by M{p'), it is easy to see that the products over A; e "D in (47) cancel 
(recall that = for A; e V) and we get finally 

4-3. Summation over p' in (21) 

In (38), (40) and (49) we have obtained all factors for the calculation of the normalized 
matrix element in such a form that the dependence on the summation indices p' is 
explicit: 

\ °l 0/ V i&v lev J 

}■£ ^ n + 6 (-i)«n + {-ly-rra ' ^ ^ 

(GI? meT>,m>l 

where the p'-independent factors IZ'^ can be read off from (38), (39) and (40), (41). The 
superscript u stands for n odd and even, respectively. 

Collecting all factors which depend on pi, I E V, the problem of performing the 
multiple summation over pi, I eV, reduces to the calculating following sum (proved in 
the Appendix): 

.tTev Ui<m,i,mevii-iy'ri + {-iy-r^) ^ ^ 

ae.(2n/a)//--)/^,r^)/^ 



X 



= < 



Ca [b ± a) (Ujev ^ ' T «^ ) — ^-ff iTh — ^ ' ^ ' 

cp{l^ ah) [±ah{{.^j, + 1 j — ^| (±h~) — ' ^' " ' ^^^^ 



with D = \V\, 

and we abbreviate 

fi = e'«'' - 1 , = e^^'' - , = 6^^''+'^- - 1, 

a = - /3 = 1 - a^ft^ (53) 
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In the calculation of the matrix element (21), we restrict ourselves to the case 
(Tq = cr^ corresponding to D odd. The case ctq 7^ (Ttt corresponding to D even can be 
done similarly. Also the two cases of even and odd parity of n have to be considered 
separately. For odd n, taking into account (38), (51) and using (28) with 

(-l)-a(a^ + l)(-C +{-irs^){{-ira + h){Y{,^^e^^^-{-irah) - 

-{-ira{a^ - 1) (C + i-irso) {{-ira - b) {Ui^^e''' + (-l)'^°a&) = (54) 



we have finally: 








.2 

n 




( ^0,0 


^0,0 ) 



kev 

^ ll;gP ^ne^ fi Qi LLl<m,l,m.eV yl ^ 'rn) 

ni<m,i,mevi-^y'hi,m H/ei^ (^ + (-l)^°«,i) ' 

where 

Analogously, in the case of even n and ctq = cr^, using (38) we get the same formula 
(55) for the matrix elements but with 

Cn={b+{-ira) a-i(a/3)(2-)/2. 
5. Product of matrix elements 

In this section we sketch the calculation of the conjugate matrix elements ($i|(J^|$o); 
where the vectors and |$o) shall have the same eigenvalues as the vectors |$i) and 
($o| used in the previous sections. This calculation can be performed in the same way 
as we did in Section 4. In analogy to (21), (22) we get for the homogeneous case 





\o 


.z 
n 




( ^0,0 


^0,0 ) 



with 



2ro ^ 



nr=i (-C + (-1)"'^) nr=i (c + {-ly^n) 

Now we have to make the same transformations as we made for R{p') in Section 4.1. 
The expression for R*{p') is obtained from R{p') given by (22) just by substituting 
b ^ —b (in particular, ( —(). Let us compare Rq{p') and Ro{p'). From the solution 
of the Baxter equations it follows that 

Q'^^'\pi) Qf^'\pi) = Q]:^'\pi) Qf^'\pi) unless lev and pi =1. 

Q'^^'\i) Qf^'\l) = - Q'^^'\i) Qf^'\l) if lev. 

So, in the final formula we have to substitute qi — > 



(^>1 


o 


z 
n 


$0) 


(*0,0 


*o,o) 
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Using these rules, in the case of ctq = cr^, from (55) we get 

roCnCa Yl{-{-lfb + ark)x 
kev 

where h^^^, fi, are h^m: fh 9i from (53) with the replacement qi R* is R 

with b — > —b. 

In the product of (55) with (56) nice simplifications appear. In R • R* we can use 

llii-l fb + ar,) {-{-!)' b + ar,) = 11:^^' (^7) 

k€V kev 
and /; ■ /; = gi ■ gl = A{qi), so that 



($0 





z 
n 




a 


z 
n 


$0) 


(^0,0 


^0,0) (^0,0 


^0,0) 



kev 

^ 4rf A(gO^-^ -|-|- (n + r-mf ^ 

in + (-i)-oa,-J (n + (-i)'^°«-,-J \hi,m\' 



(58) 



where c* is c„ with b —b and c„c* = " "is independent of parity of n. Using 
further (44) and the short notations: Aq — {—iy°So, = (— l)'^''s^ with ctq = a^^: 

i> ,2 ~ sin l{qi + q.,n) rf - rj sin |(g, + g„) 

\hi,m\ = 2 (cosg„ - cosg;) . t.^ — = — A{qm) A{qi) 



sin|(g/-gm) sini(g,-g^)' 

= -| , 2n(n + Ao)(n + A.) = (Ao + A.)(n + (-i)'^"«J(n + (-i)'^°«-,-J , 

we find that in (58) the factors A{qi) combine to 

ni(b n U^'f- n - n - f ■ 

kev 'eP V m>i,me2? ^^'V ze^u© ° 

Since also the factors a and (3 in (58) can be collected as follows {\V\ + D — n — 1): 

we get finally for arbitrary n and ctq = o-j, 





0-^ 

n 


$1) ($1 


(7^ 
n 


$0) 


(^0,0 









4i (Ao +n) (A, +n) 



X 



X n r_^_sin^iq^ 

l<n:i%ev''^ 2(?/ +1m) 

where 5=1. In a similar way we can find the product of matrix elements in the case 
of (To 7^ (^n- The final result is (60) with 5 = 0. Observe that using (60), the explicit 
appearance of excitations of type (ii), i.e. k E V has disappeared from our formula 
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(recall that we still exclude k E V). We still have normalized our matrix elements by 
the norm taken from the auxiliary system. In the next section we shall normalize to 
the norm of periodic states so that the spin matrix element becomes independent of the 
special normalization of the periodic states chosen. We shall also include the hitherto 
excluded case (iii). 

6. Final formula for the squcire of the matrix element 

The proper quantity to consider for the matrix element of the spin operator, which does 
not depend on the normalization of the eigenvectors of the transfer-matrix, is 

($o|<|^i>(^iKl^o)/((^o|^o)(^i|^i» . (61) 

The factor 

($o|$o)($i|$i)/(^o,o|^o,o)' (62) 

providing this change of normalization has been derived in formulas (74) and (75) of [2] 
for odd and even n, respectively. To obtain (61) we just divide (60) by (62). The final 
result for the square of the matrix element (60) as well as the formulas for the squares 
of norm (62) were obtained for the eigenvectors with eigenvalues not containing factors 
(A — Sq^,y, i.e. up to now we have excluded the case (iii) in Section 3.2. 

6.1. Excitations j eT> producing factors (A — Sq.y in t^{X). 

Now we explain how to modify these formulas for the matrix elements and norms if 
the eigenvalue polynomials of states |$o) and |$i) contain factors (A — Sg.y for some 
j, so that t^{rj) = for p satisfying (—1)'' = (—1)^^^. We denote the set of such j, 
corresponding to both states |$o) and |$i), by X>. As already mentioned in Section 3.2, 
case (iii), we cannot normalize the solution to the Baxter equation for the state by 
Q^'^^^\o) = 1. However, we may normalize it by Q^'^^'^\l) = 1. For the other state, 
the solution to the Baxter equation is managed by the case (i) and there is no problem 
with the normalization Q^'^^'^\o) = 1. We shall use the normalization Q^'^^'^\l) = 1 
for all j eT> irrespective of whether we have case (i) or (iii). 

In [2] the norm (56) and the matrix element (21) was calculated using the 
normalization Q^ '^'^'^O) = 1 for all j = 1, 2, . . . , n — 1. These formulas were obtained 
for the case of generic parameters for which any normalization is possible. Let us trace 
the changes in these formulas if instead we choose the normahzation Qy^^^\l) — 1 for 
j from a subset "D C {1, 2, . . . , n — 1}. 

Compare the formulas for the matrix element (21) corresponding to a different 
normalization of the solutions to the Baxter equation, i.e. we compare the term 
corresponding to a set p' in one formula with the term corresponding to the set p' in 
the other formula (the set p'^^ is obtained from the set p' by shifts pj pj + 1, j G P). 
For N = 2 this just means interchanging pj = and pj = 1, while the other components 
remain unchanged. Also, we change all rj — > — r^, j e X> in the second formula. 
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From the Baxter equations (7), (8), the solution Q^'^^^\pj + 1) normahzed by 

Q^'^^f\l) = 1 coincides with Q^'^^^'Vj) with the normahzation g^'^^''^(0) = 1. 

The only factor which is changing in (21) is N{p'). The denominator is unchanged under 
the simultaneous substitutions Vj — > —Vj and pj — > + 1- The change in the numerator 
is corrected by the division of the matrix elements corresponding to solutions of Baxter 
equations with Q^^^''\l) — 1, not by (^o,o|^o,o) but by (^^q p+g). From the 
general expression for the norm, (20) of [2] at A'^ = 2, the change of normalization 
means multiplying our matrix elements by 

(*o,o+si%+5) nr<J.(^-(-i)^-+n(-i)N' 

where = 1 if m e "D and Pm — ^ otherwise. Finally, the factor (—1)"^' gives the sign 
(-1)1^1". The formula for the norms undergoes the same changes. Therefore formally 
the expression (61) has to be invariant with respect to the substitutions rj — >■ —rj. But 
the final formula for (61) is given for the case of Ising model where we already replaced 
the dependence on Sg^ using the coincidence with rj. So the substitutions rj —rj in 
this final expression for (61) given in terms of Aq, and r^. A; = 1, 2, . . . , n — 1, will also 
change the eigenvectors entering (61) to the eigenvectors corresponding to eigenvalue 
polynomials with factors (A — r^)^, j E V, instead of (A + r^)^. This is exactly what we 
need. 

Summarizing: The final result (61) was obtained for eigenvectors with eigenvalues 
not containing factors (A — Sq.^ and which could be given in terms of Aq, A^ and rfc. 



\2 



k = 1,2, . . . ,n — 1. Now, if the eigenvalue polynomial contains the factors (A — s 
instead of (A + Sg.)"^ for some j, we just have to replace rj — > —rj in the final formula 
for all such j. 

6.2. Final result in terms of Aq, X-n, ^k, and qi. 
So the final formula for the matrix element becomes 



($0 


< 




< 




($0 


$o)(*i 


$l) 



t,i \ri-rm sm ^(gi + Qm) 



X 



A„ Ukodd,ieven(i-rk+ri){+rk-ri)) 

X —. ■ ,(64) 

2^nfe6©(+2rfc) Yl^^^^^^^^^^(^^+rk+ri){-rk-ri)^ Uk<i,k,ieyen{{+rk+ri){-rk-ri)^ ' 



where 



An = T^ — ;rr w ... — ^-t^ — ^rv^ — ^, for odd n. 



Hfeg^^i) i^o+Tk) Ilkevw (A^ - rl) U^^^w (A^+r^) UkevW iK - 4) ' 



nfcgp(°)(Ao+^fc)(A7r+rfc) 

" " (Ao + A.) n,,p(i) (Ao+r.)(A.+r,) Ukevi^^ (Ag - rl){Xl - rl) ' 
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P = P U P, = U V^^^ = U and ±r„ is the short notation for 
rm ii rn & T), ±rm if m G P and — if m G respectively. The right hand side of 
the first hne is just the right hand side of (60) except for the last product in the first 
line of (60), which is partly cancelled by the terms in the first line of the change of 
normalization (62) taking from (74) and (75) of [2]. Also, we have taken into account 
the contributions oi k e V: these have the same form as those for /c G V, just with 
reflected Vk — > — as discussed above. 

We claim that these formulas prove the formula for the matrix element which was 
given in [4] in an equivalent version. In order to show the equivalence, in the following 
we perform some transformations of (64) to get a formula which is more appropriate for 
comparison. 

6.3. Final result in terms of momenta 

Let {qi, q2, ■ ■ ■ , q/f } and {pi, p2, . . ., Pl} be the sets of the momenta of the excitations 

presenting the states |$o) from the NS-sector and from the R-sector, respectively. 
After some lengthy but straightforward transformations of (64) we obtain 









< 




(*o| 







NS R 
K pNST-r2 T\T L pBT-\2 AT y-rK i-rL 

where NS/2 (R/2) is the subset of quasi-momenta from NS (R) taking values in the 
segment < q < tt , NS/2 (R/2) containing with odd k (even A:): 



sp sin^' {sl-sl){sl-sl)' 

. NS NS „ R 



^ _ Sg + Sp _ riq^ jSQ + S^) riq^ rip^ (^q + 5p) 

^^a,P 1 J R NS 



Fornodd: .s ^ ^ £^^^ q Pf = P ^ 0, 

,R_oNs_ 1 ,_np^(^^+5p) 



Sq)(So + Sq)' ' (^tt + Sp) (sq " Sp) ' 



pit _ piNS _ ~ T _ '-'-V ^ q 

~ TT ~ „ I „ ' ~ „ NS "J 1 



+ So 



Hq- 



s. 



^J 



for n even: ^ns ^ fq pR ^ ^p p ^ tt 

" (Stt + Sq)(SO + Sq)' " (Stt " Sp) (Sq - Sp) ' 

NS 

pR _ _ pNS _ 1 T _ rip^ (^TT + gq) 

~ TT ~ ' S "-^ ■ 
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7. Bugrij— Lisovyy formula for matrix element 

In [3, 4], the matrix elements of between eigenvectors of symmetric Ising transfer- 
matrix 

tisingSym = exp (^ELi (^l-i f^fc) exp {Yl=i K ^1) exp (lELi (^k) ■ (66) 

were given. Since (25) and (66) are related by a similarity transformation with 

exp (I ELl^a;(T^-l(Jfc), 

which commutes with a^, it is natural to compare (65) with the square of the matrix 
element as given in [4]: 

I Ns( qi, ^2:■■■,^K\o■^\ Pi, P2, • • • , PL )rP = 



^ HnL 7(qfc)+7(9) L i-rR ■ . 7(Pi)+7(p) /- i^(K-L)y2 

t\ ^ Hp sinh 11 n sinh ^^HilllM ' - f-^ 



X 



^ - 2 qfc-qfc/ i ^• 2 Pi-Pj/ „;„},2 7(qfc)+7(Pi) 

n ^''' TT ^ TT ^^^^ 2 

„• u2 7(qfc)+7(qfe/) 1 1 „ • . 2 7(R)+7{Pii) 11 sin^ 3^ ' ^ ^ 

A;<A;' Sl^J^ 2 Kl' ^^^^ 2 l<k<K ''^^^ 2 

i<i<z. 

In this formula the states are labelled by the momenta of the excitations, and the squared 
matrix element is given ior a^, m — 1, . . . ,n. The operators with different values of 
m are related by similarity transformations with the translation operator. The states 
under consideration are eigenvectors of the translation operator [15] with eigenvalues 
which have unit absolute value. This explains why the formula presented does not 
depend on m. Therefore it is sufficient to calculate the matrix element for a^. The 
factors in front of the right hand side of (67) are defined by 

/ n''^ n^sinh^ 7(q)+7(P) 

C = ( (sinh 27^3; sinh 2ir J -2 - 1)^/^ = — \ „ — 5 ti-ttt 

where 7(q) is the energy of the excitation with quasi- momentum q: 

cosh 7(q) = _ " ■ cos q , (68) 

whereby (26) — tanhXj; — ab, ty — tanhi^^y = {a — b)/{a + b). 

The excitation with quasi-momentum q leads to the multiplication of the transfer 
matrix eigenvalue by e"^*^'^-' in the notations of [4], and in our notation to the 
multipUcation by ±(A — Sq)/(A -|- Sq) at A = b/a. The sign is not fixed because the 
excitations arise by pairs. Comparing (68) and (27) we get 

g7(q) = . (69) 



a Sq 



Therefore 
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which leads, in particular, to 



Our next problem is to rewrite (67) in terms of Sq. We need 



(71) 



sinh 7(0) sinh 7 (tt) / (sinh 2/^3. sinh 2/^^^)^/^ ' 



ty - t-' _ Sinh 7(0) + 7(7r) _ 1 



- 1-^ sinh 2^2^ ' 2 sinh 2^^ 



The following formulas give a correspondence between different parts of (67) and (65): 

-ix 1/2 



n -t-^^ 



sinhi(7(0)+7(7r)) V^- - 

NS 

nNS ■ 7(qfc)+7(q) I pNST-r— 

q^q^^^Smn 2 _ Sq + S^ ^qfc llq/lq,.|^^q,q* 



n nf sinh liH^i+lM sinh iM+lW T]^ 

n^^p^sinh^iHii+iM ^ .Q+s, -Pp^nLi^p -p^ 

nnr«inhl(H2)±l(s) sinhlM±lW Y{f N,,,^ 



where in the last two formulas we used 



sinh^ i(7(a) + 7(/?)) = ■ ■ sin |(a - /3) • sin i(a + /3) , (72) 

X T^x ^01 



and, in particular. 



sinh^ 1(7(0) + 7(7r)) = • ^ , (73) 

together with (71) and some trigonometric identities. The formula (72) also gives 

Sinh2 2(£0±7(^ ty-t-^ , ^ 

sin^^ tx-t-^ "'^ ^ ^ 

Formula (74) is also vahd for /3 = —a, but in this case we have to use formulas (69) and 
(70). Finally, if we take into account (73) and that K in (65) is even, L is odd, then it 
is easy to see that the formulas (65) and (67) coincide. 

8. Matrix elements for the quantum Ising chain in a transverse field 

In this section wc apply the formulas for the matrix elements obtained in Section 6 to 
the derivation of the matrix elements for the quantum Ising chain in a transverse field. 
Let us start from the L-operator (24) with a — g~^^'^ and 6 = 0: 

^ / 1 + Aaf A^-i/V| \ ,^ , 
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Expanding the transfer-matrix for the monodromy matrix (2) with this L-operator we 
have: 

2A 1 " 

t„(A) = l- _7^+..., n^-- J2{alaU, + g al) , 

^ ^ k=i 

where Ti is the Hamiltonian of the periodic quantum Ising chain in a transverse field. 
Prom (27) we get the spectrum of this Hamiltonian: 

q 

where the energies of the quasi-particle excitations are 

£(q) = (1 - 2^ cos q + g'^f'^ = ({g - if + Ag sin^ ^ j , q 7^ 0, tt , 

£(0)=5-l, £(7r)=^+l. 

In (76), the sign +/— in the front of £(q) corresponds to the absence/presence of 
the excitation with the momentum q. The NS-sector includes the states with an 
even number of excitations, the R-sector includes the states with an odd number of 
excitations. The momentum q runs over the same set as in (27). Then the formula for 
matrix elements for a^^ is given by (65) with Sq = g/e{a). After some simplification we 
get the analogue of (67) now for the quantum Ising chain: 

X 



\ Z \ \ |2 -i.)- -p-|- tj'^'-' -p-|- 

|Ns(qi,q2, • • • ,qA-|cr„|pi, P2, • • • , Pl)r| =g ^ ^Ct_[_[ — 7 — rJ_J_ 




L ne{qk) ne{pi) 



2sin£^ 



(77) 



nrnp(^(q)+^(p))^ 



n;^'(^(q) + ^(q'))^nSp'(^(p) + ^(p'))^ 



and 



e 



,(q) _ ny(^(q)+^(qO) 



np(^(q)+^(p)) 

Formally, all these formulas are correct for the paramagnetic phase where g > 1 and for 
the ferromagnetic phase where < 51 < 1. But for the case < < 1 it is natural to 
redefine the energy of zero-momentum excitation as £(0) = 1 — to be positive. Prom 
(76), this change of the sign of £(0) in the ferromagnetic phase leads to a formal change 
between absence-presence of zero-momentum excitation in the labelling of eigenstates. 
Therefore the number of the excitations in each sector (NS and R) becomes even. Direct 
calculation shows that the change of the sign of ^(O) in (77) can be absorbed to obtain 
formally the same formula (77) but with new £(0), even L (the number of the excitations 
in R-sector) and new = (1 — g'^Y^'^- 
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Formulas (67) and (77) allow to rcobtain already known formulas for the Ising 

model, e.g. the spontaneous magnetization [16, 17]. Indeed, for the quantum Ising 

chain in the ferromagnetic phase (0 < (7 < 1) and in the thermodynamic limit 

n ^ 00 (when the energies of |vac)Ns and |vac)R coincide giving the degeneration 

of the ground state), we have — 1 and therefore the spontaneous magnetization 
Ns(vacK|vac)R = C^/2 = (l-^2)i/8_ 

9. Conclusions 

In this paper we calculated the normalized spin matrix element between arbitrary states 
of the Ising model, the main result being the formulas (64), (65) and (77). We started 
with the result (12) obtained in our previous paper [2] using the Sklyanin method of 
Separation of Variables by which we obtained explicit wave functions in terms of the 
solutions of Baxter equations. The result (12) was obtained for the general N — 2 
BBS- model which is related to a generalized free-fermion Ising- type model (3). For 
this general model we were able to get the much simpler formula (21) which, however 
still involves a multiple sum over intermediate spins. Since for the general model we 
cannot perform the summation, for the further calculation we restricted ourselves to 
the Ising model with parameters (23). Performing a number of technical manipulations, 
we succeed in calculating the multiple spin sums explicitly. Although the intermediate 
formulas are quite involved, a number of surprising cancellations take place which lead 
to the rather simple formula (60) for the spin matrix element square. This comes still 
normalized to the auxihary states involved in the method of Separation of Variables, but 
it is not difficult to convert this result into the properly normalized matrix elements for 
the model with periodic boundary condition. The final formula becomes more lengthy 
due to normalization factors. We show by which transformations we get the recently 
conjectured formula of A. Bugrij and O. Lisovyy. Our derivation provides a first proof of 
these formulas. Another application of the formulas obtained in this article is the result 
(77) for the spin matrix elements for the finite quantum Ising chain in a transverse field. 

The presence of degenerations in the spectrum for the special Ising parameter values 
forced us in this case to normalize the Baxter equation solutions differently for different 
excited states. The complexity of the formulas gives httle hope that for more general 
parameter values the multiple spin summations can be done in the near future, even if 
then the degeneration problems may be avoided. 
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A. Proof of the summation formula 

The aim of this Appendix is to find a factorized expression for the sum 

where is given by (48) with (46). The result for = ±0;^^ is 

{-^r {{-^Tn ± ((-i)'"^ + C) 



/ X n r2n/a) f(^-i)/2 {D-z)i2 

^ ^ lll,meV,l<m \P^'''l;m) ^^^^ 

cir° (a ± 6) (a n,,. e.^^ T n,e. (^n/a) (/^«,)°/'-' ^ ^^^^ 



= < 



Y\.l,meV,l<m i^^l,' 

with D^\V\, 



.odd 
-13 



and we use the abbreviations (53). Replacing ±aq^ — > ±Pq^ in (79), see (46), (48), 
amounts to ol ^ (3 in c°^'^ and cg|^^°, and 

(6 ± a) ( n,g7, e'^^ Tab) ^ (1 T a 6) ( 6 n^ep e'^^ ± «) ^ odd, 

(a ± b) (a rijei, e'^'^ ^ 6) ^ (1 ^ a 6) (±a 6 [1^62? + 1) D even. 

In Section (4.3) the last bracket in the numerator of (78) is also needed with 
+(^ = b/a replaced by —C,. However, in order not to complicate the formulas, we shall 
always use (78) as it is written here. Since the only other 6-dependence, which is in 
and Tj, is quadratic, this can be adjusted at the end of the calculation just by changing 
the sign of b. Also, it is sufficient to prove formula (79) for the case — aq^. Then the 
results for the cases = — cc^j and = ^Pgi can be obtained by simple transformations 
of the variables r; and b. 

First, we find the recurrence relation for Yj) with respect to D = It relates this 
quantity for D and D—2. Then we verify (79) for small D. Finally, by explicitly inserting 
our conjectured solutions (79) into the recursion relation we prove the expression. The 
cases of D odd and even have to be treated separately, since we are dealing with a 
two-step relation. 

A.l. Derivation of the recursion relation 

This recurrence relation for Yx> is obtained by using the identity, compare (17), 
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when the number of Xi exceeds the number of Uk at least by two. Fix any index s E T> 
and consider D + 1 values of Xi and two values of yf {yi, 2/2} — {—Cs, ~C}) 

{xo, X,,..., xd} = {ts, -r„ -i-iyWi, -i-irr2, . . . , -{-1^^ , -{-lY'^rD}. 

omitted 

Since two parameters yk are chosen, we must have D > 3. Now we separate the two 
terms in (81) which correspond to i = 0, 1 and present them as a summation over 
Ps G {0, 1} for {xq,Xi} = rs(— 1)^^ Then (81) becomes 
(-l)'"(r,(-l)'"+?.)(r,(-l)"- + 



E 



Ps 



^ (rl - rl) a^,,,(-r,(-l)^^ + n(-l)N ' ^ ^ 

Now in (78) we separate the summation over a certain fixed discrete variable Ps and 
use the identity (82) to replace the summation over ps G {0, 1} by a summation over k. 
After this, we move the summation over k to the front of formula (78) and collect the 
factors depending on p^. 

Yii^s^ki-^ririi-iy +ii){ri{-ir + Q 



^ , -iy''(rk(-iy>'+^k)(rk{-iy''+ 2rsi-rk{-iy'^+Q{-rk{-iy'^+ c) 

X 



n,^,.,M-i)" + ri(-l)«) {rl - r?) n,^i,,(-r^(-l)« + r,(-l)") ' 

Now after multiplication in the second line of (83), we perform the summation over 
Pk e {0, 1} by means of the relation 

E,, i-iyini-iy^ + e,)(-rfe(-l)^'= + Q = -2rk{^k - Q ■ 

Finally, we get the desired recursion relation, valid for any fixed index s & T>: 

V _ \^ V ^^fc^^ - rl) (6 - ^s) ,x 

l^^g \'k ' s) \.\.li^k,s\' I ' k) 

where Yx)/k,s is just the original Yx> with the indices k and s removed from T). 

Since (84) is a two step difference equation, in general it has two independent 
solutions Y.^'^ and Y^^ . The quantity (78) is equal to a linear combination of these 
solutions, Yx) = ciY^^ + 02 Y^^ and the constants ci, C2 are fixed by calculating Yx> 
explicitly for small D. 

In the following, using (84), we give the details of the proof of (79) for the case D 
odd. The derivation for D even proceeds analogously. 

A. 2. Proving the summation formula for D odd. 

For D = 1 we start calculating directly (78), = Og^: 

Yis} - E(-l)"(-^(-l)" + -s)(rs{-ir + hi a) = 2rs^^^±0:^^ , 

ps ^ ' 
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which proves (79) for this particular case. This result, together with the recursion 
relation (84), defines Yjy for odd D uniquely. For example, for the set V — {si, S2, S3}, 
the recursion relation (84) with selected s = S3 gives 

^ 4rir3 {C -rl){ai - Q3) 4r2r3 (C^ - r|)(tt2 - 0:3) 

^W,.2,.3} ^{^2}^^2_^2) (rl-rf) + ^{^i> (r2 - rf) (r? - r^) 

(6 + a) ULi C^rifi) fe'^^ie"^^ - ab)h2,s ^ e'^^{e'^' - ab)h,,3 



a' (-/3) ni<R„^<3 hi,m V (e^^'^ - (e'^"^ - e'^') 

_ + «) n?=i (2n/0 ,^3 



- aM-^)n../... (nt.e--a.), (85) 
where we used 

ak-as^ -a pk,s, ( - ^ , ri -r^^ pi,k 86) 

Ik 9k Ji Ik 

with pk^s — (e'^* — e}^'')/{gk9s)- Observe that the big bracket in the second fine of (85) 

factorizes and leads to a result symmetrical in the three indices. The result obtained in 

(85) proves (79) for D = 3 and ^; = aq^. We can easily continue this recursive procedure 

to conjecture the formula (79) for odd D. To prove it, it is enough to show that the 
right-hand side of (79) satisfies the recursion relation (84). The right-hand side of (79) 
can be presented as Yj) = ciY^^ + C2Y.^\ where ci = —{a -\- b)ab, C2 = a-\-b and 



n/n I \ AD-l)/2 (£>-3)/2 
=c^ , YS^=Y^'^Y[e^^^. (87) 

1 \.l<m '^l,m I 



Thus to prove (79), it suffices to prove that Y^^^ and Y^^^ satisfy (84). 

Let us prove that Y^'' satisfies the recurrence relation (84) with = aq^. Using 
the elementary relation c'^^2/'^1d'^ — —ot^~^/3^~'^ and (86), we reduce the problem to 
the proof of the identity 

Tf / \ AD-l)/2 {D-Z)/2 TT /9„ /_\ /.(D-3)/2 (D-5)/2 

\{i{2ri/a)fl g\ ^ ^ lhjLk,s{^ri/a) fi g^ ^ 



X 



^TkVs fsft^gt'e'^" rr figi 



a 



2 



n 



l^k,S 

We see that all vi and the fi, gi, hi^rn with /, m ^ k, s match between both sides of (88), 
leaving us with 

which is equivalent to the interpolation identity (81) used here for the following choice 
of the parameters: Xk = e'^* for e P/s, = e~'^', 

{2/0, 2/1, • • • , Vd-z] = {0, a^,...,a^ , a'^, . . . , a'^} . 




The proof that Yj^ satisfies the recurrence relation (84) reduces to the same 
interpolation identity, but with uq omitted. This proves (79) for D odd and = aq^. 
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As was explained before, the change of the sign at C, in (78) can be adjusted by 
changing the sign of h in the final formula. Similarly, the change of the sign at 
can be adjusted by the simultaneous change of the the signs of ri and Finally, the 
transformation h 1/6 leads to ri ri/b'^, — ttg, Pqjb'^- This allows to find (78) 
^-t ii = Pqi if we know (78) at = —cxqi- The mentioned transformations cover all the 
cases of (78) for D odd. 
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